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ABSTRACT 

Context. Local-box hydrodynamical model atmospheres provide statistical information about a star's emergent radiation field which 
allows one to predict the level of its granulation-related micro- variability. Space-based photometry is now sufficiently accurate to test 
model predictions. 

Aims. We aim to model the photometric granulation background of HD 49933 as well as the Sun, and compare the predictions to the 
measurements obtained by the CoRoT and SOHO satellite missions. 

Methods. We construct hydrodynamical model atmospheres representing HD 49933 and the Sun, and use a previously developed 
scaling technique to obtain the observable disk-integrated brightness fluctuations. We further performed exploratory magneto- 
hydrodynamical simulations to gauge the impact of small scale magnetic fields on the synthetic light-curves. 

Results. We find that the granulation-related brightness fluctuations depend on metallicity. We obtain a satisfactory correspondence 
between prediction and observation for the Sun, validating our approach. For HD 49933, we arrive at a significant over-estimation 
by a factor of two to three in total power. Locally generated magnetic fields are unlikely to be responsible, otherwise existing fields 
would need to be rather strong to sufficiently suppress the granulation signal. Presently suggested updates on the fundamental stellar 
parameters do not improve the correspondence; however, an ad-hoc increase of the HD 49933 surface gravity by about 0.2 dex would 
eliminate most of the discrepancy. 

Conclusions. We diagnose a puzzling discrepancy between the predicted and observed granulation background in HD 49933, with 
only rather ad-hoc ideas for remedies at hand. 
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1. Introduction 

The bright (my = 5.77) F-type dwarf HD 49933 (HR2530, 
ADS 5505 A) was the first prime seismic target of the CoRoT 
space photometry mission (" Auvergne et al.l2009h . A clear signa- 
ture of solar-like oscillations wa s detected in th e acquired light- 
curve dAppourchaux et ^ 120081: iMTchel et al.l l2008). This con- 
firms expectations about the presence of solar-like oscillations in 
HD 49933 which were anticipate d from results o f ground-based 
observations of radial velocity (Moss er et alj|2005l) . 

It is now generally accepted that solar-like oscillations are 
excited by stochastic convective motions taking place in the stel- 
lar surface layers on time-scales of the same order as the periods 
of the excited modes. The convective flows are thermally driven, 
and imprint a time-varying brightness pattern - granulation - on 
the stellar surface. While the oscillations manifest themselves as 
discrete peaks in a temporal power spectrum of a star's light- 



* The CoRoT space mission, launched on December 27th 2006, has 
been developed and is operated by CNES, with the contribution of 
Austria, Belgium, Brazil, ESA, Germany and Spain. 



curve, the stochastic evolution of the granulation pattern leads to 
a continuous signal, the so-called granulation background, since 
it usually constitutes the background on which oscillatory peaks 
are superimposed in a spectrum. Oscillations, as well as gran- 
ulation, induce very small brightness fluctuations, challenging 
the sensitivity level of photometry achievable from the ground, 
and motivating photometry from space such as that performed 
by CoRoT. 

The modelling of stellar granulation - mainly driven by 
the need for improved stellar atmosphere models of late-type 
stars - made substantial progress over the last two decades. It 
has now reached a level that should allow a direct comparison 
between the predicted and observed granulation properties, in 
particular the granulation backgrou nd (Trampedach et al. 1998t 
ISvensson & Ludwigll2005l:lLudwigll2006.: Guenther et al. ,20081) . 
Hitherto, the comparison was somewhat hampered by the limited 
photometric accuracy, duty cycle, and duration of existing stellar 
measurements. The outstanding accuracy achieved on HD 49933 
by CoRoT offers a new opportunity for a stringent comparison 
between predicted and observed granulation background, and is 
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attempted here. We will see, however, that some problems are 
encountered. 

The paper is organized as follows: we briefly describe the 
photometric data (Sect. |2]i, and the hydrodynamical granulation 
models (Sect.O, comment on how we derive the synthetic back- 
ground spectra (Sect. 13. It , present comparisons for the Sun and 
HD 49933 (Sect. |4]i, and finish by discussing effects which may 
play a role in the encountered problems (Sects. |5]and|6]l. 



2. CoRoT photometric data and power spectra 



We us ed reduced photometric data of the N2-level jSamadi et al.l 
12007 ') covering 137 d for HD 49933 from long run LRaOl, 
and 157d for the F-dwarfs HD 181420 and HD 181906 from 
long run LRcOl. W e applied the same bolometric coiTection 
dMichel et al.ll2009l) to all targets, amounting to a factor of 1.11 
in amplitude of the brightness fluctuations. We calculated power 
spectra by applying a standard Fast Fourier Transform from the 
bolometric light-curve without further processing. In particular, 
we did not apply any coiTection to the long-term trend apparent 
in the light-curve of HD 49933. The reason was that applying 
the measured loss of instrumental sensitivity due to aging again 
produced a trend but in the opposite sense. This in turn resulted 
in a power spectrum almost indistinguishable from the power 
spectrum of the uncorrected data. On the other hand, simply de- 
trending by removing any linear long-term variation had an im- 
pact on the power - but at frequencies v < 0.3 mHz only. These 
frequencies are not relevant to the conclusions drawn in this pa- 
per. Hence, we did not further investigate this issue, but remark 
that there is an instrument-related uncertainty in the measured 
power and its slope in the range v < 0.3 mHz. 

When we refer to "power" in this paper we always mean 
power spectral density. We follow the normalization that the 
power integrated over the interval between zero (excluded) and 
the Nyquist frequency equals the variance cr^ of the signal in the 
time domain. 



3. CO^BOLD radiation-hydrodynamics simulations 

We used the radiation-hydrodynamics code CO'' BOLD (for 
further informatio n about the code and applications see 
iFrevtag et al.l 120021: IWedemever et al.l 1200 4) to construct three 
3D Cartesian "local-box" model atmospheres. Two models were 
intended to represent HD 49933, and one - for reference - 
the Sun. Table [1] summarizes the model properties. The atmo- 
spheric parameters for HD 49933 were chosen close to val- 
ues reco mmended in the literature , T^ff - 678 + 130K 
dBruntt et al . 2008), and logg * 4.25 dAppourchaux et al. .2008.). 
The metallicity of the star is [Fe/H] « -0.37 dex dSolano et al.l 



120051: iGillon & Magainll2006l) . The opacities used in the models 
for HD 49933 are based on detailed low-tem perature opacities 
from the MARCS stellar atmosphere package dGustafsson et al.l 
I2OO81) . We only had these data available for particular metal- 
licities, and not for arbitrary element mixtures. Due to this re- 
striction, we decided to compute two 3D models with metal- 
licities which bracket the star's observed metallicity. As will 
be demonstrated below, the chemical composition has an im- 
portant impact on the model properties. The solar model Sunl 
is a CO^BOLD model of an older generation e mploying opaci - 
ties from the ATLAS model atmosphere package (Kuruczl l2005h . 
While this leads to some systematic differences between the 
models, we believe that they so not influence the comparisons 
performed in this paper. Our two hydrodynamical model atmo- 



spheres for HD 49933 wer e also employed to predic t excitation 
rates of oscillatory modes (Sama di et al.ll2008ll2009l) . 

Unfortunately, the aforementi oned st ellar parameters are not 
completely reliable. iBruntt et akl d2004l) list effective tempera- 
tures and surface gravities from various determinations which 
differ by almost 500 K and 0.5 d ex, and the debate about the 
exact parameters is still ongoing (lBruntd|200^; 'KaflingeretiD 
2008). It may turn out that HD 49933 is 200-300 K cooler than 
assumed in our models. More important for the level of the 
granulation background, however, is the the value of the sur- 
face gravity. We will come back to this point later but want to 
emphasize that HD 49933 has an accurately known parallax of 
33.7 + 0.4mas from HIPPARCOS ( vanLeeuwen 2007) which 
helps to constrain its surface gravity and stellar radius. In this 
paper, we use a radius of 1.35 fo r HD 49933, within the er- 
ror box given by Appo urchaux et all ([2008) of 1 .36 ± 0.06 Rq. 

Since the brightness fluctuations which are central to the 
present investigation are intimately linked to the heating and 
cooling provided by the radiation field to the gas, we give here 
some details about the treatment of the radiative transfer in the 
hydrodynamical models. All models employ 140x 140x 150 grid 
points for the spatial mesh in the x-, y-, and (vertical) z-direction, 
respectively. For all models, the mesh spacing is equidistant in 
the horizontal directions. In vertical direction, SO and SI em- 
ploy a non-equidistant spacing, Sunl again an equidistant mesh. 
In the important layers around t = 1 the spacing in vertical di- 
rection amounts to 46.9 km (0.19 Hp"'', definition see Table [TJ, 
48.2 km (0.20 i/^"'), and 15. 1 km (0. 1 1 H'p'"^) for models SO, S 1, 
and Sunl, respectively. Note that the models SO and SI use the 
same vertical mesh. The difference in resolution at t = 1 is 
the result of how the model adjusts on the computational mesh. 
The radiative transfer is solved along long cha racteristics em- 
ploying Feautrier's method (e.g.. iMihalasll 1 9781) assuming strict 
LTE. Corresponding "rays" are started at each point of the spatial 
mesh in the top-most layer. All models use bundles of rays in one 
vertical and two inclined directions, as well as four azimuthal di- 
rections (M - 3 case described in Ludwig 2006) coinciding with 
the coordinate axes. In total this amounts to 176400 rays along 
which the radiative transfer equation is solved at each time step. 
The number of equivalent wavelength points Nobm to represent 
the wavelength dependence of the radiation field is five for mod- 
els SO, and Sunl, six for SI, as given in Table [T] The necessary 
interpolation of the radiative heating or cooling between the sys- 
tem of rays and spatial mesh (for hydrodynamics) is performed 
in an energy-conserving fashion. To optimize performance, the 
treatment of the radiative transfer is switched to the diffusion 
approximation in the deep, optically thick layers. 

Finally, we want to point out two aspects which are of par- 
ticular importance for the interpretation of the granular back- 
ground. Qualitatively, the lower metallicity in model S 1 in com- 
parison to SO leads to a lower overall opacity (primarily due to 
the lower H" opacity at lower electron pressure). This in turn 
implies an increase of the mass density at given optical depth. 
At optical depth unity, we find an increase of the mass density 
by 40 % in model S 1 relative to SO. At fixed flow geometry, 
the higher density would require a smaller temperature contrast 
or smaller velocity differences between up- and down-flowing 
material to transport the prescribed stellar energy flux. What ac- 
tually happens in SI is an increase of the temperature contrast 
reflected by an increase of the intensity contrast (see Table [Til. 
This is likely a consequence of the dependence of the opacity 
on temperature. On the other hand, as expected, the typical ve- 
locity amplitudes in model SI are smaller by about a: 10% at 
the maximum of the vertical velocity. Finally, the typical gran- 
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Table 1. CO^BOLD radiation-hydrodynamics model atmospheres: "Model" is the model's name used in this paper, Teff the effective 
temperature, log^ the gravitational acceleration, [Fe/H] the metallicity, R an assumed stellar radius (not intrinsic to the simulation 
proper), /x = ly the linear horizontal size of the square-shaped computational box, /z its vertical extent, T the duration of the simulated 
time series, Tc the sound crossing time over //p"^, H^'^ the pressure scale height at Rosseland optical depth unity, the 
relative spatial white light intensity contrast at disk-center, o-gi-^n the temporal relative disk-integrated granulation-related brightness 
fluctuations in white light (in parenthesis an estimate of the uncertainty), Nobm the number of equivalent frequency points considered 
in the solution of the radiative transfer equation, "Modelcode" an internal identifier of the model sequence. 
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SO 


6720 


4.25 


0.0 


1.35 


16.4 


24.2 


135.6 


29.5 


0.241 


0.196 


84 (17) 
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d3t68g43mm00n01 


Si 


6730 


4.25 


-1.0 


1.35 


16.1 


24.2 


73.8 


29.1 


0.238 


0.229 


64(13) 
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Sunl 
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4.44 


0.0 


1.00 


5.6 


2.5 


52.2 


17.8 


0.141 


0.172 


40(8) 


5 


d3gt57g44n53 



ular scales become shorter in S 1 . The combination of these fac- 
tors determines the amplitude and characteristic frequency of the 
granulation background. 

3.1. Synthetic power spectra of the observable flux 

The hydrodynamical models provide time series of the radi- 
ation intensity at different limb -angles includ ing the temporal 
convection-related fluctuations. Ludwigl (l2006h showed that, to- 
gether with an assumption about the stellar radius, this infor- 
mation is sufficient to predict the power spectrum of the disk- 
integrated, observable, relative brightness variations. He ob- 
tained for the frequency component of the power spectrum the 
relation 



(/r) 
iff 

where 



{^m=l ^fnf^m (^m)j 



(1) 



(2) 



Angular brackets denote expectation values, / the observable 
disk-integrated flux, and / its (complex) Fourier transform. The 
asterisk indicates the conjugate complex, c„, a discrete integra- 
tion weight, jj.^ the cosine of the limb-angle, and /,„ the spatial 
average of the intensity in direction cosine ju,,,. According Eq. ([T]l 
the observable power scales inversely proportional to the num- 
ber of simulation boxes of surface area A tiling the visible 
stellar hemisphere of radius R. For a detailed derivation and dis- 
cussion of Eq. O see lLudwig|(l2006l) . 

Due to the finite duration of a simulation, the obtained power 
spectra of the relative brightness fluctuations exhibit appreciable 
noise. To improve the effective signal-to-noise level and facil- 
itate the comparison to observations we fit a simple analytical 
model to the simulated spectra. The analytical model reflects 
our prejudice that the frequency dependence of the granulation- 
related background signal should show rather little structure. 
IVIoreover, it allows us to eliminate the acoustic eigenmodes of 
the computational domain which are excited by the convection- 
related fluctuations but are not directly comparable to the ob- 
served modes. We applied a model for the s pectral power d en- 
sity P which is a generalized Harvey model (lHarvevlll985h for 



the background, plus a sum of Lorentzians for the box-modes 
according to 



P(y) = 



ak 



(3) 



V is the cyclic (temporal) frequency. All other variables are fit- 
ting parameters of the model: b is the asymptotic amplitude of 
the background towards low frequencies, vi and V2 characteris- 
tic frequencies of the background, ci and a2 power law expo- 
nents, are the mode amplitudes of K considered modes, 
their frequency positions, wa their line widths. We performed a 
maximum likelihood estimation of the parameters; we used the 
commonly adopted model that each frequency component of a 
power spectrum is statistically independent and follows a prob- 
ability distribution p which is a -distribution with two degrees 
of freedom. This is an exponential of the form 



p(x) = -^exp(--^l. 

(x) \ {x) 



(4) 



{x) is the expectation value of the random variable x. The 
model (O allows us to extrapolate the fitted spectrum to fre- 
quencies higher than the Nyquist frequency VNyquist- To mitigate 
effects of aliasing present at the highest frequencies of the simu- 
lated spectra, we did not fit P{v) itself but P(v) + P(2vNyquist - v). 
This includes the "mirroring" of high frequency power into the 
frequency domain of interest by aliasing. As we will see, P{v) 
decreases strongly towards high frequencies so that aliasing ef- 
fects are only noticeable close to the Nyquist frequency. 

To maximize the likelihood function we used the POWELL 
function implemented in IDL. The maximization problem turned 
out not to be well conditioned, and some manual intervention 
was always necessary to obtain a stable and acceptable fit. This 
points to the need for a more sophisticated fitting procedure with 
suitable regularization. However, for the present purposes we 
considered the fitting rather a kind of constrained smoothing so 
that the shortcomings from the viewpoint of statistical stringency 
were not considered vital. We verified that our fitting resulted in 
a representation close to what was obtained by direct smoothing 
of the simulated spectra. 

4. Results 

4.1. Fitting the simulated raw spectra 

Figure [T] illustrates the result of fitting the analytical model of 
Eq. (O to the simulated raw spectra calculated from the time 



4 



Ludwig et al.: The enigmatic granulation background of tlie CoRoT target HD 49933 



Table 2. Fitted parameters of Eq. ^ for the synthetic power 
spectra of the three hydrodynamical models. Background- and 
line-amplitudes are given in ppm^/y^Hz, frequencies in mHz. 
Missing values for a line indicate that it was left out since un- 
necessary to obtain a satisfactory fit. 



Parameter 


Hydrodynamical models 




SO 


SI 


Sunl 


b 


6.056 


2.881 


1.046 


Vl 


0.9593 


1.259 


1.295 


ai 


2.486 


3.434 


2.953 




2.383 


2.821 


3.667 


012 


8.061 


9.278 


8.035 


ai 


0.4130 


0.5389 


2.149 


«i 


1.266 


1.674 


3.252 


Wi 


0.09507 


0.08464 


0.02861 


ai 


4.498 


1.732 




U2 


1.676 


2.027 




W2 


0.01155 


0.02994 




fl3 


2.188 


2.142 




"3 


2.107 


2.618 






0.01217 


0.004683 





series with the help of Eq. ([T]). Table |2] lists the fitted parame- 
ters. As already mentioned, the simulated spectra show an ap- 
preciable level of noise. In the figure, the simulated spectra were 
smoothed to allow a clearer comparison. By comparing solar 
models with different numerical set-ups, we estimated that the 
absolute uncertainty in the simulated spectra amounts to about 
+20 % in total (i.e. frequency integrated) amplitude; we expect 
the relative accuracy among models of similar numerical set- 
ups to be somewhat better Note that we left out the contribution 
of the box modes in the fitting functions in Fig. [T] The simu- 
lated spectra show clear differences in the total power between 
the Sun and HD 49933 - primarily reflecting the change of the 
granular cell size relative to the stellar radius among the stars. 
The metal-depleted model for HD 49933 SI shows only 58% of 
the model SO in total power This is the result of a reduction 
of the granular scale at lower metallicity, over-compensating of 
the increase in the contrast (cf. Eq. ||5l). Moreover, the power 
in model SI is shifted towards lower frequencies. This we trace 
back to longer granular life-times in SI, presumably related to 
the lower typical velocities in SI. All models predict a rather 
steep drop of the granular background signal towards high fre- 
quencies, similar to the one observed in the Sun. 

4.2. Comparison with VIRGO observations of the Sun 

Figure [2l co nstitutes an update of results presented by 
ISvensson & L udwig (2005), and shows a comparison (not fit) 
between the disk-integrated, photometric fluctuations derived 
from the solar model Sunl, and observational data from the 
VIRGO instrument on board the SOHO satellite. Our focus is 
the high frequency region of the solar signal in which granular 
contributions dominate. The VIRGO power spectrum was cal- 
culated from (level 2) time series data provided by the VIRGO 
team. One year of data (1996.5-1997.5) was chosen close to so- 
lar minimum activity to minimize the possible contribution of 
activity related variability. The plotted power spectrum is based 
on data of the green channel of the VIRGO three-channel sun- 
photomet er (SPM, see|F r6hlich et al. ( 1997)). Following the pro- 
cedure of ISvensson & Ludwig. (,2005) . it has been converted to 



lO.OOp 




0.01 0.10 1.00 10.00 

Frequency v [mHz] 



Fig. 1. Inter-comparison of the raw (solid lines) and fitted 
(dashed lines) spectra (power spectral density) of the three hy- 
drodynamical simulations. The contribution of box-modes to the 
fitted signal was left out in the graphs. 



1 0.000 F 




0.001 I ." , 

0.01 0.10 1.00 10.00 

Frequency v [mHz] 



Fig. 2. Power spectra of disk-integrated photometric fluc- 
tuations for the Sun: the predicted background signal of 
model Sunl (green/grey solid line) and observational data from 
SOHO/VIRGO (black solid line). Note the steep decline in 
power towards high frequencies. 



white light fluctuations by matching (by shifting in power) a 
corresponding power spectrum based on VIRGO PM06V-A ab- 
solute radiometer data in the frequency range 0.3 to 2.0 mHz. 
Instead of using the PM06V-A power spectrum directly, this 
rather involved procedure was necessary since PM06V-A and 
SPM power spectra deviate substantially in the high frequency 
region. The authors found little information about this mismatch 
in the literature. However, it appears to be agreed that the SPM 
spectrum reflects the actual solar behavior dFrohlich et alJll997t 
lAndersen et alJ[T998h . in particular showing the steep (roughly 
as v"'*) decline at the highest frequencies. 

We find a satisfactory agreement of the continuous back- 
ground signal between model predictions and observations be- 
tween 0.1 and 8.0 mHz. In particular, the background in the 
p-mode frequency range is matched quite well. The observed 
background power shows a "kink" at frequencies slightly above 
1 mHz relative to the model prediction, and the "plateau" power 
around 0.2 mHz is not exactly represented. iMichel et al.l (l2009h 
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interpreted the "kink" as the sign of the transition between a 
meso-granular and granular contribution to the power spectrum. 
This interpretation is obviously not supported by our granulation 
model, which already provides too much power around 1 mHz, 
and has a geometrical dimension which would not allow the de- 
velopment of meso-gra nular convection cells of 5- 10 Mm size 
("Novem ber et al.iri981h . Moreover, Michel and co-workers at- 
tribute a similar total fluctuation amplitude to meso-granulation 
and granulation. According to Eq. (|5]l, this would correspond to 
a spatial brightness contrast which would be a sizable fraction 
of the granulation contrast and should be easily recognizable in 
images of the solar surface - which is not the case. 

Towards low frequencies (v < 0.1 mHz), the increase of 
the observed signal is due to the solar magnetic activity, and 
cannot be represented by the pure (non-MHD) hydrodynamical 
model applied here. We also considered the power spectrum at 
solar maximum (year 2000.5-2001.5, not shown). Towards so- 
lar maximum, the activity -related power increases but leaves the 
granulation-related signal intact; the power spectrum is essen- 
tially identical at frequencies v > 0.2 mHz. It is difficult to pro- 
vide an objective criterion to determine to which lowest tempo- 
ral frequency our model can make reliable predictions. A limit is 
set by the maximum size that features can develop in the simula- 
tion. However, it is not obvious how this relates to their life-time, 
hence the frequencies on which they will have an impact. 

A conspicuous feature exhibited by the observations as well 
as the model is the "step-like" shape of the granular background 
signal. As it will be shown later, this feature is not obvious in the 
observations of HD 49933, while it is predicted by the simula- 
tions. 



4.3. Comparison with CoRoT observations of HD 49933 

Figure |3] shows a comparison between the predicted power spec- 
tra (power spectral density) of models SO and SI. Ideally, the 
observed spectrum should fall somewhere between models SO 
and SI. This is not at all the case; the predictions lie noticeably 
higher with a power almost a factor of two greater than observed 
around 0.9 mHz. Moreover, the overall shape of the observed 
spectrum is not well represented by the models. In particular, 
the conspicuous drop of the background level is not obvious in 
the observations. 

In Fig. [3] we added two spectra of F-dwarfs which were also 
observed by CoRoT. The stars have parameters not too different 
from HD 49933 (HD 181420: Feff ^6650 K , log ^^4. 17 [ M/H] ^ - 
0.04; HD 181906: 6380/4. 15/-0. 14; see Michel et al.l (|2008|) ). 
The overall levels of the brightness fluctuations difi'er, however, 
the spectra show a remarkable similarity in shape. We take this 
as an indication that HD 49933 is a typical representative of its 
class, and not an exceptional case among F-dwarfs. We note 
here that the similarity of the spectral shape at y > 0.3 mHz 
might be connected to the ra ther similar rotational peri ods of the 
stars of 3.4 d for HD 4993 3 (Appourc haux et al.ll2008l) . 2.8 d for 
HD 181906 (.Garcia et al.ll2009,) ^ and 3.5. . .4.5 d for HD 181420 
dBarbanet al.1l2009l). 



Michel et al.l (l2008l) assumed that the continuum in the 



p-mode region indeed corresponds to the granulation back- 
ground, fitted it with a standard Harvey model, and arrived at 
granulation-related brightness fluctuations of CTgjan ~ 40ppm 
with appreciable uncertainty. This would put the fluctuations 
in HD 49933 on the same level as observed in the Sun, and 
is significantly lower than our predictions. However, accept- 
ing for the moment a much lower level of fluctuations, we ex- 
perimented in an ad-hoc fashion by adding (besides photomet- 




'VlJJiIi, |HD181906 




0.10 1.00 
Frequency v [mHz] 

Fig. 3. Power spectra for HD 49933: dashed lines depict the pre- 
dictions from models SO and SI. The solid lines are again the 
predicted spectra but increased by a constant corresponding to 
the observed photometric white noise level. The black line de- 
picts the spectrum of HD 49933. The other two solid lines depict 
spectra of the similar F-dwarfs HD 181420 and HD 181906. 




0.10 

Frequency v [mHz] 



Fig. 4. Like Fig. [3] but with scaled models SO and SI (for de- 
tails see text) and adding an ad-hoc signal for magnetic activity 
and photometric noise (solid lines). The dashed lines depict the 
models after scaling without added ad-hoc components. 



ric noise) an activity-like signal (assuming a power-law) to our 
simulated granulation spectra, and simultaneously scaled their 
power and characteristic frequency to enforce correspondence 
with the observations. Figure|4]illustrates the outcome after scal- 
ing the power derived from the models by a factor of 0.33 and 
the frequency by a factor of 1.6. This resulted in a plausible - 
albeit not entirely satisfactory - coiTespondence. Later we shall 
try to interpret the scaling in terms of a change of fundamental 
stellar parameters (see Sect. 15. lb . 



5. Discussion 

Can we find a reason for the mismatch between model and ob- 
servation for HD 49933? In the following we discuss the impact 
of uncertainties in the fundamental parameters and the influence 
of magnetic activity on the mismatch. 



6 



Ludwig et al.: The enigmatic granulation background of the CoRoT target HD 49933 



5. 1 . Uncertainties in the fundamental parameters 

As alluded to before, the fundamental parameters of HD 49933 
(in particular Teff, logg, and R) might be noticeably different 
from the values we assumed in our modeling. Ludwig (2006) 
give an approximate relation between the spatial contrast of the 
granulation and the photocentric variability, which can be equiv- 
alently stated in terms of the relative temporal (RMS) brightness 
fluctuations cr„y^„ as 



R 



(5) 



/gran IS the typical linear size of a granular cell, R the stellar ra- 
dius, / the (spatially resolved) emergent intensity, 5/,ms its spa- 
tial standard deviation, and I its spatial average. We wrote a 
proportionality only since we want to investigate the differen- 
tial functiona l depen dence on the fundamental parameters here. 
iFrevtag et al.l (11997 *) argue on the basis of numerical simulations 
that /gran is proportional to the pressure scale height at the stellai^ 



surface //^"'"^ 



'gran ^ '^p 



(6) 



When writing the last proportionality, we assumed a constant 
mean molecular weight. For the total (frequency integrated) 
power o"gj.aj, of the brightness fluctuations, it follows that using 
the fundamental relationship between luminosity L, radius, and 
effective temperature 



2 eff 
C^gran 



g2R2 



eff 



5hn 



(7) 



Introducing the luminosity in the second proportionality in 
Eq. (|2l) was motivated by the notion that the accurately known 
parallax and visual magnitude of HD 49933 essentially fixes its 
luminosity. Relation (|7]) is not exact since the scaling of the gran- 
ular size with stellar parameters is not exactly accounted for; 
for instance, the formula does not accurately represent the scal- 
ing of CTgran among the models as given in Table [1] However, 
we think it is good enough to provide an estimate of the dif- 
ferential effects in vicinity of the stellar parameters we used in 
our modeling. As evident from Table[T] the granulation contrast 
does not sensitively depend on atmospheric parameters among 
F- and G-dwarfs, making it largely invariant over the discussed 
interval of temperatures, gravities, and chemical compositions 
for HD 49933. Keeping the luminosity fixed gives estimates of 
the change o f crl^,„ by -hO.58 dex w ith the recently suggested pa- 
rameters o fKallingeret all ( l2008l) (r,ff^645 K, log ^=3.9), and 
-1-0.048 dex with the ones of iBrunttI (|2009|) ("evolutionary so- 
lution" reff=6560K, logg=4.19). Despite the high exponent of 
Teff in Eq. (|7]i, the changes are primarily driven by the changes 
in log g. We are left with the situation that the suggested param- 
eters would increase the discrepancy between models and obser- 
vations when applied in the modeling. For the set of parameters 
suggested by Kallinger and co-workers the deterioration would 
be rather drastic. 

The ad-hoc scaling in power and frequency described in 
Sec. 14.31 to enforce the correspondence to observations was 
motivated by Eq. |7l and by the notion that the granular life- 
time scales with the atmospheric acoustic cut-off frequency 
(Kjeldsen & Bedding 1995; Svensson & Ludwig 2005) as v oc 
gl VT^eff- This allows us to associate the scaling factors in power 
and characteristic frequency with changes of -200 K in T^n and 



-1-0.2 dex in logg of HD 49933 (assuming constant luminosity). 
Again, while helping to reduce the present discrepancies, in par- 
ticular the increase of the gravity is at odds with the latest deter- 
minations of the star's fundamental parameters. 

5.2. Influence of magnetic activity 

By analogy to the rise of power towards low frequencies due 
to magnetic activity in the Sun (Fig. |2]l, one is tempted to at- 
tribute the rise of power seen in HD 49933 to magnetic ac- 
tivity as well. This is quite plausible, since the lig ht-curve 
of H D 49933 exhibits clear signatures of spottedness ( Mossed 
2009). Moreover, HD 49933 is a rather rapid rotator, favoring dy- 
namo action. However, polarimetric observations (with NARVAL 
and ESPADONS) indicate that there is no magnetic field present 
above a: 1 G whi ch is structured on spatial scales of the order of 
the stellar radius dCatala & al.ll2009l) . This raises the question of 
what kind of mechanism could produce the apparently ubiqui- 
tous field on smaller scales? 

There is mounting theoretical and observation al evidence 

dCatt aneo 1999; Trujiflo Bueno et al. 2004; Vogler & Schiissled 
2007J) that dynamo action is possible in convective surface flows 
generating magnetic fields on very small scales. To investigate 
whether such magnetic fields can explain the puzzling shape 
of the spectrum around 1 mHz in HD 49933, we performed 
explo ratory MHD simulations with CO^BOLD ( Steffen et al.l 
l2009h . They were restricted to two spatial dimensions, and per- 
formed for solar atmospheric conditions. We calculated several 
runs with different, prescribed magnetic flux levels correspond- 
ing to field-free to weak plage regions. Figure |5] shows tempo- 
ral power spectra of the light curves of the simulation runs (not 
scaled to full disk). Despite that the simulations were tailored 
to solar conditions, we think that the very similar topology of 
granular flows among solar-type stars makes them qualitatively 
also applicable to HD 49933. The power spectra of the light- 
curves show a trend towards lower power at higher frequencies 
for increasing total magnetic flux. Obviously, small scale mag- 
netic fields alter the convective dynamics on small spatial scales 
but do not lead to the formation of structures on large scales 
which may leave an imprint at lower frequencies in the bright- 
ness fluctuations. The results are preliminary, and further work 
is necessary. However, for the time being we consider it unlikely 
that local dynamo action in the granular flow - if at all present 
— can explain the shape of the power spectrum of HD 49933 in 
the p-mode frequency region. 

Equation (|5]l shows that a reduction of the granular scale 
at fixed stellar radius could reduce the observed granulation- 
related brightness fluctuations. It is known that magnetic fields 
of sufficient strength have this effect (for a recent example see 
IJacoutot et al.l l20()8i and references therein). However, to ob- 
tain a sizable reduction one would need a mean field strength of 
several 100 G. A local dynamo appears unlikely to sustain such 
magnetic flux levels (Vogl er & SchiissleJ 12007 ). which leaves 
us with the need for an efficient generation of magnetic fields on 
larger scales. 

6. Conclusions 

3D hydrodynamical model atmospheres can be applied to predict 
granulation-related disk-integrated brightness fluctuations. We 
demonstrated that the theoretical approach successfully repro- 
duces the convection-related background signal in the frequency 
region of the observed solar p-modes. Unexpectedly, we found 
a mismatch between predicted and observed background when 
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Fig. 5. Temporal power spectra of the horizontally averaged 
emergent intensity in vertical direction of different 2D MHD 
runs assuming different levels of initial field strength (solid lines 
of different colors/shades of grey). The simulations were per- 
formed assuming solar atmospheric parameters. 



applying the same approach to HD 49933, which we argued is a 
normal representative among the F-dwarfs observed by CoRoT. 
We excluded local dynamo action in the granular flow as the cul- 
prit for the mismatch. The existence of a magnetic field for other 
reasons with a large filling factor and sufficient strength (several 
lOOG) able to noticeably influence the granular dynamics could 
qualitatively explain the observed low level of granular bright- 
ness fluctuations. Whether this is viable from the viewpoint of 
dynamo theory is unclear. In this context we discussed the ob- 
servational finding that the star, while magnetically rather active, 
does not exhibit a m agnetic fi eld organized on scales comparable 
to the stellar radius (ICatala & al. 2009). 

Recently sugges ted new sets of fundamental stellar param- 
eters for HD 49933 (iKalUnger et"aLl|2008; Brunt3l2009l) are un- 
likely to remove the discrepancy - primarily since they suggest a 
lower surface gravity than we assumed in the present modelling. 
In contrast, when increasing the gravity by 0.2 dex (to about 
logg=4.45), and decreasing the effective temperature by 200 K 
(to about 6550 K) we can enforce a plausible, albeit not per- 
fect, correspondence between theory and observation. Whether 
this is compatible with constraints from stellar structure, 3D 
atmospheric modelling, and p-mode excitation needs to be as- 
sessed. For the moment we are left with a puzzling discrepancy 
between the predicted and observed granulation background in 
HD 49933, with only a rather ad-hoc idea for a remedy at hand. 

The granulation-related brightness fluctuations are sensitive 
to the stellar metallicity. A reduction of the metallicity from solar 
to I/IO solar reduced the total fluctuation power by a factor of 
0.58 at the studied Teff and log g. 

If in the future one is able to devise an observational or the- 
oretical calibration of the granular background signal, its high 
sensitivity to the surface gravity could render it an effective 
gravimeter. CoRoT has obtained high-precision photometry for 
many giants. The stronger convection-related background sig- 
nal makes this easy to detect and quantify. We consider giants 
as additional, interesting test cases for our understanding of the 
dynamics of convective flows on the surface of stars. 

Acknowledgements. The authors thank Michel Auvergne for providing infor- 
ination on the change of instrumental sensitivity, and Hans Bruntt for coininu- 
nicating his latest results on the fundamental parameters of HD 49933 before 



publication. HGL would like to thank Claude Catala for discussions about spec- 
troscopy and polarimetry of HD 49933 on several occasions. HGL acknowledges 
financial support from EU contract MEXT-CT-2004-014265 (CIFIST) 



References 

Andersen, B., Appourchaux, T., Crommelynck, D., et al. 1998, in lAU sym- 
posium 181: Sounding Solar and Stellar Interiors, ed. J. Provost & F.-X. 
Schmider (Kluwer Academic Publishers), 147 
Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 
Auvergne, P., Bodin, P., Boisnard, L., et al. 2009, ArXiv e-prints 
Barban, C, Deheuvels, S., Baudin, P., et al. 2009, submitted to A&A (this vol- 
ume) 

Bruntt, H. 2009, submitted to A&A (this volume) 

Bruntt, H., Bikmaev, I. P., Catala, C, et al. 2004, A&A, 425, 683 

Bruntt, H., De Cat, P, & Aerts, C. 2008, A&A, 478, 487 

Catala, C. & al. 2009, in preparation 

Cattaneo, E 1999, ApJ, 515, L39 

Ereytag, B., Holweger, H., Steffen, M., & Ludwig, H.-G. 1997, in Science with 
the VLT Interferometer, ed. E. Paresce (Springer), 316-317 

Ereytag, B., Steffen, M., & Dorch, B. 2002, AN, 323, 213 

Erdhlich, C, Andersen, B. N., Appourchaux, T., et al. 1997, Sol. Phys., 170, I 

Garcia, R., Regulo, C, Samadi, R. Ballot, J., et al. 2009, submitted to A&A (this 
volume) 

Gillon, M. & Magain, P 2006, A&A, 448, 341 

Guenther, D. B., Kalhnger, T., Gruberbauer, M., et al. 2008, ApJ, 687, 1448 
Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 
Harvey, J. 1985, in ESA Special Publication, Vol. 235, Future Missions in Solar, 

Heliospheric & Space Plasma Physics, ed. E. Rolfe & B. Battrick, 199 — 1- 
Jacoutot, L., Kosovichev, A. G., Wray, A., & Mansour, N. N. 2008, ApJ, 684, 

L51 

Kallinger, T., Gruberbauer, M., Guenther, D. B., Fossati, L., & Weiss, W. W. 

2008, submitted to A&A 
Kjeldsen, H. & Bedding, T. 1995, A&A, 293, 87 
Kurucz, R. L. 2005, MSAIt Suppl., 8, 14 
Ludwig, H.-G. 2006, A&A, 445, 661 

Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558 
Michel, E., Samadi, R., Baudin, E, et al. 2009, A&A in press 
Mihalas, D. 1978, Stellar Atmospheres (Freeman and Company) 
Mosser, B. 2009, priv. comm. 

Mosser, B., Bouchy, E, Catala, C, et al. 2005, A&A, 431, L13 
November, L. J., Toomre, J., Gebbie, K. B., & Simon, G. W. 1981, ApJ, 245, 
L123 

Samadi, R., Belkacem, K., Goupil, M.-J., Ludwig, H.-G., & Dupret, M.-A. 2008, 

Communications in Asteroseisinology, 157, 130 
Samadi, R., Fialho, E, Costa, J. E. S., et al. 2007, ArXiv Astrophysics e-prints 
Samadi, R., Ludwig, H.-G., Belkacem, K., Goupil, M.-J., & Dupret, M.-A. 2009, 

submitted to A&A (this voluine) 
Solano, E.. Catala, C, Garrido, R., et al. 2005, AJ, 129, 547 
Steffen, M., Ludwig, H. G., & Steiner, O. 2009, ArXiv e-prints 
Svensson, E. & Ludwig, H.-G. 2005, in Proceedings of the 13th Cambridge 

Workshop on Cool Stars. Stellar Systems, and the Sun (ESA SP-560), ed. 

E Eavata, G. Hussain, & B. Battrick (ESA Publication.s Devision), 979-984 
Trampedach, R., Christensen-Dalsgaard, J., Nordlund, A., & Stein, R. 1998, 

in The First MONS Workshop: Science with a Small Space Telescope, ed. 

H. Kjeldsen & T. Bedding (Aarhus Universitet), 59 
Trujillo Bueno, J., Shchukina, N., & Asensio Ramos, A. 2004, Nature, 430, 326 
van Leeuwen, E 2007, A&A, 474, 653 
Vogler, A. & Schussler, M. 2007, A&A, 465, L43 

Wedemeyer, S., Ereytag, B., Steffen, M., Ludwig, H.-G., & Holweger H. 2004, 
A&A, 414, 1121 



